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We propose the reconstruction of F(R) gravity in such a way that corresponding theory admits 
cyclic and ekpyrotic universe solutions. The number of explicit examples of such F(R) model 
is found. The comparison with the reconstructed scalar-tensor theory is made. We also present 
F(R) gravity which provides the little rip evolution and gives the realistic gravitational alternative 

" for ACDM cosmology. The time for little rip dissolution of bound structures in such theory is 

estimated. We demonstrate that transformed little rip F(R) solution becomes qualitatively different 

. cosmological solution with Big Bang type singularity in Einstein frame. 
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INTRODUCTION 



It is widely accepted nowadays that universe evolution has passed through two (super) accelerating epochs: early- 
time inflation and late-time cosmic acceleration. Despite the close similarity between these two accelerating epochs 
there is a number of differences. For instance, inflationary era is characterized by large value of curvature while dark 
era occurs at small curvature. Moreover, the effective equation of state parameter w is qualitatively different for 
O I inflation or late-time cosmic acceleration. Inflation occurs right after Big Bang/Big Crunch singularity while number 
of quintessence/phantom dark energy models end up at finite time future singularity. This suggests the very natural 
conjecture: Big Bang should be identified with finite time singularity. In other words, the universe realizes cyclic 
evolution which may be produced by multi-fluid [l| or by the theory with oscillating equation of state parameter Q • 
fSJ ■ It is natural to construct the cyclic universe within the theory which describes the inflation and dark energy in unified 
^ | way. 

Modified gravity is realistic alternative for unified description of inflation with dark energy (for review of such 
unification in modified gravity, see Q). Moreover, its weak-field limit may lead to newtonian regime (for review, 
see (H). Hence, modified gravity represents very natural candidate where cyclic cosmology may be realized. The 
purpose of this work is precisely this one: reconstruction of modified F(R) gravity which leads to cyclic evolution. 
In the next section the reconstruction technique for scalar-tensor theory and F(R) gravity is developed. We find the 
explicit examples of scalar-tensor theory and F(R) gravity which have the solutions corresponding to cyclic universe. 
In third section F(R) action is presented as General Relativity plus ideal fluid. The reconstruction is developed for 
such formulation of modified gravity. Again, the example of F(R) gravity which admits cyclic universe solution is 
found. In fourth section using above methods we show that ekpyrotic scenario may be also realized in frames of 
modified gravity. In fifth section we reconstruct F(R) gravity which induces little rip cosmology. Such non-singular 
dark energy model was recently proposed in Ref. [5| as an alternative to ACDM model. Furthermore, we demonstrate 
that reconstruction of F(R) theory which leads to non-singular phantom evolution is possible. Some summary is given 
in the last section. 

II. CONFORMALLY CYCLIC UNIVERSE 

Let us show how to construct the cyclic universe solution in modified gravity. The starting point is the spatially 
flat FRW metric 

ds 2 = -dt 2 +a{t) 2 (dx'f . (1) 

i=l,2,3 
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where a(t) is the monotonically expanding function of the cosmological time t. Let a(t) has the following form 

lno(t) = H t + h(t). (2) 
Here Hq is a positive constant and h(t) is a function with period T: 

h(t + T) = h(t). (3) 

Then 

ds 2 = -dt 2 +a(t + T) 2 ^ ( da 

1=1,2,3 

= -dt 2 +a(t) 2 e 2H ° T J2 ( d ^~ 



!x") 2 



^ix 

i=l,2,3 
, 2 



= -dt 2 +a(t) 2 ( di T ■ ( 4 ) 

i=l,2,3 

Here x 1 = e H ° T x l . Eqs. (JTJ) and (01 show that the metric at t is physically identical with that at t + T. In this 
sense, the universe is cyclic although the universe is monotonically expanding. Hence, it called as conformally cyclic 
universe. A simplest example is 

h(t)=H lC os(^j . (5) 

Here Hi is a positive constant less than Hq, which guarantees that the universe is monotonically expanding. 

We now consider the models which reproduce the scale factor ait) given by ([2]) with ([5]). First, let us investigate 
the scalar-tensor model, whose action is given by 

S = J d 4 X^(^R - ^{(f)d^d^ - V(4>) + Matter j . (6) 

Here, w((/>) and V(<p) are functions of the scalar <p. For the model where oj(4>) and V(</>) are given by a single function 
f{4>), as follows, 

«W = -4/"(0) , V(0) = 4 (3/'(0) 2 + , (7) 

the exact solution of the FRW equations has the following form: 

= t, H=f'(t). (8) 

In case of fl2J with (jSj, one finds 

/(*) = In a(t) , f(t) = H - ^sin . (9) 

Therefore if we consider the model 

2Hi /2tt\ 2 /2tt<A 



u ™ = —{t) co \t j 



the conformal cyclic universe ([5]) with ([5]) can be reproduced. 
We now consider F(R) theory 



•V, = j d A Xyf=g + ^matter j ■ I 1 I i 
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Here F(R) is an appropriate function of the scalar curvature R. The action (fTTj) is equivalently rewritten as 



S = J + W))+Anatter} 



(12) 



Here, P and Q are proper functions of the auxiliary scalar <j>. By the variation over </>, it follows that = P'((f>)R+Q' (</>), 
which may be solved with respect to <f> as <f> = <f>(R). By substituting the obtained expression of MR) into (fT2")l . one 
arrives again at the F(R) gravity action (for reconstruction of modified gravities, see Refs. 0, 

S = J d^xj—g | ffl^ + £ raattor J , F(R) = P((f>(R))R + Q(<t>(R)) . (13) 

If the scale factor a is given by a proper function g(t) as a — aoe 9 ^ with a constant ao, by solving the following 
second-rank differential equation, 



= 2 ^7^ - MW^T + ^"Wm + E (1 + «*) W^e- 311 ^* , (14) 



we can find the form of P{4>) which gives a solution a = a^e 9 ^. One also finds the form of Q{4>) as follows: 
Q(<t>) = -6 (g'Wf P&) - Gg'W^f + E P>oao 3il+W * ] e-^ +w ^ . 

i 

In case of ((23) with |5j), when we neglect the contribution from matter, Eq. (fT4l) has the following form: 



Since one cannot solve (|15[) explicitly, we assume P(</>) can be expanded by the Fourier series: 

P(0) = f>„cos(^) + f: ^sinf^V (16) 



n=l 



Then one gets (ft = pi = and for n > 2, 

2 (rc - l) 2 4ttP (n - 1) n - 4 

P " ~ ~ff 1 (n + 2) Pn - 1_ T> + 2) tl + ^2 P "" 2 ' 

2(»-l) 2 47rg (n-l) n-4 

ill (n + 2) 2 [n + 2) n + 2 

In p7|) , we put q — 0. Eq. IT7|) shows the existence of the solution of (fT5l) and the corresponding F(R) gravity model. 
When matter can be neglected, Eq. (fT4"]) can be rewritten as 

i(/(»r"") = >-^, da) 

which gives 

In the second equality, we have used the partial integration. Furthermore by writing P(<f>) as 

P(0) = C/(0)~ 4 , (20) 

(flT)]) is rewritten as follows: 



2 dU{4>) 12 [,,{dU(4>)\ 



(21) 
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As an example, we may consider 

U(<j>) = Ui + U 2 cosu)<t>. (22) 
Here Ui, J7 2 , and ui are constants and it is assumed C7i > U 2 > 0, which show U{4>) does not vanish. Hence 

,,,x _ 2U 2 uj sinc^ 3L/|^ 2 (20 - 1 sin (2gg) + C) 
5 C/i + [/ 2 cos^ ([/j + [/ 2 cosw^) 2 

Here C is a constant of the integration. Note that the obtained g'(4>) is not always positive. We also note that the 
g'(4>) is not periodic due to the term proportional to cf>. This kind of term always appears for the periodic U((j>) since 
the integrand of the last term in (|2ip is positive definite and therefore the integration cannot be periodic. 
Let us now consider an example where the function P((j>) is given by, 

P(<f>) = U((j))- 4 = P o (cosw0) 4 , (24) 

where Po and ui are constants. Then, by the equation (|18j) . the solution is, 

g'{4>) = ffo(cosw0) 2 + 2w(shi2w(j) — tanw</>) , (25) 

where go is an integration constant. Note that here the solution contains divergences due to the term of the tangent 
in (|24p . These divergences correspond to points where the scale factor becomes null aito) — 0, which can be identified 
with a Big Bang/ Crunch singularity. This class of singularities is very common in cyclic Universes, where the ekpyrotic 
scenario is reproduced. In the same way other cyclic universes may be used for explicit reconstruction of modified 
gravity. 



III. RECONSTRUCTED F(R) AS AN EFFECTIVE PERFECT FLUID 

Let us start with the modified gravity FRW equations written in the following form: 



3H 2 = 



1 



F'(R) \2 



-F(R) + 3Hd t F'(R) -3H 



3H Z -2H = - 



1 



F'{R) \2 



-F(R) + 2Hd t F'{R) + d 2 t F'{R) - H 



(26) 



where we have neglected the contributions of any other kind of matter. If we compare Eqs. (|26|) with the standard FRW 
equations (3H 2 — n 2 p and —3H 2 — 2H — n 2 p), we may identify both right sides of Eqs. (f2l)|) with the energy-density 
and pressure of a perfect fluid, in such a way that they are: 



F'(R) \2 



F(R) + 3Hd t F'(R) -3H 



P = 2 



' ~F(R) + 2Hd t F'(R) + d 2 F'{R) ) + H 



(27) 



_F'(R) \2 

Then, Eqs. (|2^|) take the form of the usual FRW equations, where the EoS parameter for this dark fluid is defined by 



F'(R) 



'}F(R) + 2Hd t F'{R) + d 2 F'{R)) + H 



■pfa (±F(R) + 3Hd t F'(R)) - 3H 



The corresponding EoS may be written as follows: 



1 



AH 



' " 2 -'^^ H d t F'(R)\ . 



F'(R) 



d?F'(R) 



F'(R) 



(28) 



(29) 



The Ricci scalar is R — 6(2H 2 + H), then F(R) is a function of the Hubble parameter H and its derivative H . The 
form of the EoS is written as: p = —p + g(H, H, H, ■ ■ ■ ), where 



g(H, H, H, ■ 



4H + d? {In F'(R)) + {d t \nF'(R)) 2 - HdtlnF'(R) 



(30) 
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Then, by combining the FRW equations, it yields the following differential equation: 

H+^-g(H,H,H,---) = 0. (31) 

Hence, for a given cosmological model, the function g given in (|30p may be seen as a function of cosmic time t, and 
then by the time-dependence of the Ricci scalar, the function g is rewritten in terms of R. Finally, the function F(R) 
is recovered by the expression (|30|) . In this sense, Eq. (|3Tj) combining with the expression (|30|) results in 



dt 



+ x(t) 2 - H(t)x(t) = 2H(t) , (32) 



where x(t) — '"^^ ■ In general, this equation is very difficult to be solved for a particular Hubble parameter, 
so that the reconstruction of F(R) gravity is not so simple. Let us consider an example to illustrate this method, we 
assume a power-law evolution for the Hubble parameter, 

Hit) = j . (33) 

Then, the solution of the equation (|33[) is given by, 



„ -1 + Vl-6a + a 2 (-1+ . 2Cl 2 



Recalling that x(t) — d ( ln F J^ R ^) ^ ^ e F(R) action yields, 



F(R) = dR J fi? ^ +c 2 \ , (35) 

where {Ci,C2} are integration constants. Hence, the F(R) action has been reconstructed for the class of power-law 
solutions expressed in Q33p. However, for a particular oscillating solution, as the ones studied in the section above, 
the equation (|32|) turns out very complicated, and in general, one can not obtain an explicit expression for the action. 
Nevertheless, let us proceed in a different way, by reconstructing the Hubble parameter assuming a particular x(t) in 
the equation (p?2"j) . 

, . d(kiF'(R(t)) , . 

x(t) = — ^ v " = Hi cos cot . (36) 

Then, by the equation (|32[) . the Hubble parameter obtained is given by, 

H(t) = Cie"" 1 ^""* +e" ffl 2-" t / die" 1 *""* (Hf + HI 2 coscot - 2H!UJsmujt) . (37) 



F(R) = I dR e-ir^M(ii)] _ 



And the action is obtained, 



Here the function t(R) comes from the inversion of the definition of the Ricci scalar R(t) = 6(2H 2 + H). Hence, this 
kind of cyclic evolution (|3"T|) can be reproduced in the frame of F(R) gravities. However, the complexity of the integral 
does not allow us to get an explicit expression for the F(R) action. Note that for explicit values of the constant Hi, 
and of the frequency u>, some information may be obtained for different times, since each term in the integral (|37[) 
becomes important for different values of the time coordinate, i.e. at different epochs of the Universe evolution. 



IV. THE EKPYROTIC SCENARIO IN F(R) GRAVITY 

The inflationary epoch (early-time acceleration) just after the initial singularity that was the origin of the Universe, 
is able to solve some of intrinsic problems of the Big Bang model. These problems, known as the flatness problem, the 
question on the homogeneity at large scales, the absence of monopoles or the origin of the inhomogeneities that lead 
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to the formation of large scale structure, can be solved quite well by the inflationary paradigm. However, a decade 
ago a new proposal, alternative to inflation, was suggested, the so-called Ekpyrotic/cyclic Universe (see Ref. [9|). This 
new scenario solves additionally the puzzles of the standard cosmological model as well as it provides perhaps a more 
complete picture of the Universe evolution. While in the inflationary scenario, the initial conditions are required to 
complete the cosmological picture, as the time begins at the Big Bang singularity, in the ekpyrotic Universe, this is 
not a requirement due to its cyclic nature. The evolution presented by this class of cyclic Universe contains in general 
four stages per cycle: a first initial hot state similar to the standard Big Bang model, then a phase of accelerated 
expansion, after which a phase where the Universe starts to contract occurs and finally ends in a Big Bang/Crunch 
transition, when the cycle starts again. The corresponding period where the main problems enumerated above are 
solved occurs during the contracting phase. In the usual ekpyrotic models, a scalar field is included to reproduce a 
cyclic Universe. However, it is clear that modified gravity, and precisely F(R) gravity, can perfectly reproduce the 
ekpyrotic scenario. In order to show that the above models, reproduced by F(R) gravity, can solve the initial problems 
by means of a cyclic Universe containing a contracting phase, let us consider a FRW Universe described by, 

3 rr2 PrnO PrO PaO fc /„„, 

—H = —3- + — + — g 2 + Pf(b) ■ (39) 

K a° a* a? a 

Here p m Q is the energy-density for pressureless matter, p r Q for radiation, p a o for the anisotropics, k is the spatial 
curvature of the Universe and Pf(r) is the effective energy-density defined in the previous section for the extra 
geometrical terms. Then, in order to get a homogeneous and isotropic spatially flat Universe in a contracting phase, 
the effective EoS parameter WptR) > 1, so that when the scale factor tends to zero, the F(R) terms in the equations 
dominate over the rest, and the result is the same as in the inflationary scenario. 

Looking at the above examples, we have that for the model given in (l25l , the effective EoS parameter for the perfect 
fluid defined with the extra terms of F(R) is given by, 

Pf(r) 4w (— 2w(coswi) + oj(secu}t) 2 + Hq(cos uit) (sin u>t)) 

wf(R) = = — H 1 5 • (40) 

Pf(r) 3(£To(coscji) 2 + 2w(sin2wt)- (tanwt)) 

Hence, for each epoch of the Universe, the effective EoS behaves in a different way, we can summarize the phases of 
each cycle as follows 

W F (R) = < -|, t^TT/2uj (41) 

Hence, we can see that the F(R) fluid dominates over the rest of components in (|4TJ)) when its EoS is larger than 
1, which occurs for small times and a contracting universe, while for large times, the EoS behaves as dark energy, 
and even it crosses the phantom barrier when WptR) < — 1. So that the ekpyrotic scenario takes place, and the 
result is the observable Universe. Hence, we have shown that cyclic Universe can be well reconstructed in frames of 
modified gravity, and the initial problems, as the flatness or horizon problems, can be solved. However, generally 
the cyclic Universes contain singularities of the type of Big Bang/Crunch, where the scale factor goes to zero, and 
the energy densities for matter/radiation grow to infinity, as ,for example, occurs for the model (|2~4"1) . In the original 
ekpyrotic model, based on a scalar field, the divergences are controlled by a function introduced in the action as a 
strong coupling between the scalar field and the other components. In the kind of ekpyrotic scenario proposed here, 
a smooth transition along the Big Bang/ Crunch can be reproduced by means of the reconstruction of models that 
avoid the singularity (as it is realized for several examples given above) or by introducing a coupling in the matter 
action, 

S M = J dx A ^p{R)C m . (42) 

Hence, by an appropriate choice of the function /3(R), the matter divergences can be avoided, and the transition 
across the Big Bang/Crunch can be realized smoothly. 

Therefore, the ekpyrotic scenario can be realized with no need to introduce an additional field but only in terms of 
modified gravity. 



V. LITTLE RIP COSMOLOGY IN MODIFIED GRAVITY 



In the previous sections, several oscillating Universes have been reconstructed via F(R) gravity. Cyclic Universes 
seem to be well reproduced in these theories and the ekpyrotic scenario can be realized. Here we are interesting in the 



7 



study of cosmological models which are able to reproduce super-accelerated phase in the context of modified gravity. 
In this case, the corresponding effective EoS parameter ([251) crosses the phantom barrier, that is Wpm) < — 1. It is 
well-known that cosmological phantom models usually contain the so-called big rip singularity In the case of 

the big rip singularity, the phantom energy-density and the scale factor diverges in a finite future. One of the most 
surprising consequences of the big rip is the dissolution of bounded systems, as the Solar System and atoms before 
the singularity (see Ref . (ToL [IH ) . 

In this section, we are interesting to reconstruct F(R) gravity which is able to reproduce a phantom behavior, being 
free of future singularity. Nevertheless, such singularity free or little rip cosmology (see Refs. @, [H|]) also leads to 
dissolution of bound structures. Let us recall the first FRW equation in F(R), 

^H 2 (t) = p F{R) (t), (43) 

where the energy-density is defined in (|27[) . By definition, a big rip singularity occurs in a finite time, usually denoted 
by t s , when the scale factor a(t) and the energy-density diverges. So in order to avoid a big rip singularity, PF(R){t) 
must remain positive and finite for all t. Looking at the expression for pF(R)(t) m P7|) . it seems complicated to obtain 
a general condition on the form of F(R), but a natural condition (not sufficient) seems F'(R) > for all R in order 
to avoid divergences. On the other hand, in the absence of matter, the EoS parameter ([25)) can be rewritten as, 

2F[ 

W F(R) = -1 - ^ • ( 44 ) 

Hence, for a super-accelerated expansion, the Hubble parameter is required to be an increasing function of time. 
Looking at some of the models reconstructed in the previous sections, we can see that the Hubble parameter ([25)) 
is a periodic function that reproduces a cyclic Universe. From the effective EoS parameter (|41|). it is clear that this 
model exhibits periods of super-accelerating expansion and it does not contain any future singularity. To find out 
if the model ([25")) may lead to a little rip, one has to determine the duration of each cycle and the strength of the 
effective repulsive force reproduced by the F(R) terms and compare with the binding forces of coupled systems (as 
for example the Solar System). However, to ensure that a little rip is reproduced, one has to study models of eternal 
acceleration but free of future singularities, which may have a stronger growth in time than those models containing 
big rip singularities. Let us use the reconstruction technique using Eqs. (I12ti21[) . As an example, we consider the 
function, 

U{4>) = , (45) 

where a and ft are constants. Then, by the expression (|21[) . the Hubble parameter and the scale factor yield, 

H{t) = h e at +hi, -> a(t) = a e 4/3oQt+6Qt . (46) 



where ho = 4a(3 and hi — 6a. It is straightforward to see that the function ([46)) describes a Universe, where for small 
times t <C a, the Hubble parameter can be approximated as a constant, reproducing a de Sitter solution, as in the 
case of ACDM model. For large times, the Universe ends in an eternal phantom phase, where the EoS parameter 
w f(r) < — lj but without big rip singularity. Nevertheless, a little rip (dissolution of bound structures) might occur 
in a finite time, similarly to the model presented in Ref. [5|, as it is pointed out below. The functions P{4>) and Q{4>) 
can be easily reconstructed by the expressions ([T5j) and ([20)) . 

P(0) = e 4f3c °* , Q{cj)) = -6a 2 (3 + 4/3e^)(3 + 8/3e^)e 4/3o °*. (47) 

The F(R) action that reproduces the solution (|4"6")) can be calculated by inverting the expression of the Ricci scalar 
R = Q(2H 2 + H) and by using the expression (fT3")) . which yields, 



F(R) 



Ho 



3 y i2R " 



(48) 



where R = a 2 , d = -24e" 39 / 12 i? , and C 2 = 2 v / 3i? - Hence, we have reconstructed F(R) model which is able to 
reproduce a super-accelerated (eternal) phase, which does not lead to a future singularity. Note that the action ([4"5|) 
turns out to be the Einstein-Hilbert action plus some corrections for small values of the Ricci curvature R, where the 
exponential functions are expanded in power series, 

R 2 R 3 

F(R)~ ^R + kz— + k 315 2 +••• • (49) 
Ko K Q 



s 



Here the couplings m are constants depending on C\.i- Hence, for small values of the Ricci scalar the action reduces 
to the action for General Relativity plus power-law curvature corrections, as i? 2 , which is known recipe to cure the 
singularities and which has a viable behavior Hence, the action (|48| represents a viable model where GR can 

be recovered while the curvature scalar corrections remain small. It is important that such additional corrections 
become relevant close to the little rip evolution. In order to estimate the time for the little rip induced dissolution of 
bound structures in a naive way, one might compare the energy-density of a bound system as the Solar System with 
the density Pf(r)- For the model (|46l) . such density can be approximated for large times by 

PF(R) = Poe 2at , (50) 

where po is a constant that can be set by imposing that the current value of the energy-density is Pf(r) {to) = Jy-^o ~ 
1CP 47 GeV 4 , where the age of the Universe is taken to be to ~ 13.73 Gyrs, according to Ref. One can set the time 
of the little rip dissolution occurrence when the gravitational coupling of the Sun- Earth system is broken due to the 
cosmological expansion. By assuming a mean density of the Sun-Earth system given by p©-e = 0.594 x 10 -3 kg/m 3 ~ 
10" 21 GeV 4 , the time for the little rip dissolution of bound structures is, 

29 93 

^lr = 13.73 Gyrs H ■ — . (51) 

a 

Hence, depending on the parameter a, the appearance of the little rip may last shorter or longer. For example, when 
a = 10 _1 Gyrs , the little rip occurs at the Universe age of £lr, = 313.03 Gyrs, while for a > 1 Gyrs , the time for 
the decoupling will be much shorter. 

A lot of models can be reconstructed in frames of F(R) gravity which are able to reproduce a super- accelerating 
phase free of singularities. Let us consider, for example, a Hubble parameter that reproduces a phantom Universe 
without big rip [7j, 

H 2 = ^-(t - t Q ) 2 = H a N + H x . (52) 

Here we have introduced the number of e-foldings N = In — instead of the cosmological time t. For the Universe 
described by the function (|52|). the effective EoS parameter is less than — 1, so it also describes a phantom evolution 
but free of future singularities. By using the reconstruction technique from Ref. 0, the action that reproduces the 
above solution is given by, 

where K(a, &, x) is the Kummers serie. In this model, the energy-density is given by, 

Pf(r)= Po{t-t Q ) 2 , (54) 

where to — as the origin of the Universe evolution. Then, by adjusting the value for po with the present value of 
the energy-density as in the above model, the little rip occurs when the Universe age is t = 137.3 x 10 12 Gyrs, which 
clearly grows much slower than for above model. 

Hence, we have shown here that F(R) gravity is able to reproduce successfully a phantom scenario free of future 
singularity, where the little rip occurs at the time which depends completely on the expansion growth rate of the 
model. 

We are now interested to explore the corresponding picture in the Einstein frame for those solutions describing 
a little rip in F(R) gravity. In order to reconstruct the action in the Einstein frame, one has to use a conformal 
transformation that removes the strong coupling in the action (I12[) . 

9Ep.v = tf'gp.v , where ft 2 = P{4>) , (55) 
where the subscript e stands for Einstein frame. A quintessence-like action results in the Einstein frame 



Se = J d 4 xy/-g E 



R E --cj(4>)d^cf>-U( ( l>) 



(56) 



where 



(57) 
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while the cosmological solutions in the Jordan frame are transformed through out the conformal transformation as, 

a E (f E ) = P{<j){t)) 1/2 a{t) where dt E = P((j){t)) 1/2 dt . (58) 



As an example, one can consider the action (|47|) and (|48|) . which was found to reproduce the solution (|46|). By 
applying the conformal transformation (|55[) . the action in the Einstein frame (|57[) is described by the kinetic term 
and the scalar potential, 



w(0) = 4a 2 /3 2 e 2a0 , and [/(</>) 
The solution is transformed as, 



-6a 2 (3 + 4pe a,p )(3 + 8/3e a *)e 



0x -4/3e° 



t E = 



2/3c° 



c'- 

Z 



(59) 



(60) 



Here the second expression is given by an exponential integral. In order to calculate an exact expression for a E (i E ), 
let us approximate the exponential e at ~ 1 + at + 0(t 2 ) for small times, so that the solution can be expressed as, 



«e(*e) = OE0* E e' 



6 6/3 7 *e 



#e(*e) = 6/3 7 + 

t E 



(61) 



where 7 = ae~ 2 ' 9 . Then, it is straightforward to see that the Universe described in the Einstein frame not only 
does not describe a super-accelerating expansion evolution but it also contains a singularity, a type of initial or Big 
Bang singularity. On the same time the solution in the Jordan frame was free of singularities. The singularity of the 
Einstein frame universe is probably the manifestation of little rip (dissolution of bound structures). 

Let us consider now the inverse reconstruction by studying first the little rip scenario with quintessence/phantom 
fields described by the action ©, and then reconstruct the F(R) action by means of a conformal transformation. It is 
well-known that any cosmological solution can be reconstructed with minimally coupled quintessence/phantom fields 
scalar fields (see Ref. [15|). As an example, let us consider the kinetic term and the scalar potential 



uj{(f>) = -^hoae^ , V(cj)) = \ (3h 2 e 2a4> + (6h hi + h a)e a ^ + 3hf) . 

Using the FRW equations and Eqs.([7]), the solution for the Hubble parameter and the scalar field are 

H E (t) = h e atE + h x , <Hi E ) = ts , 



(62) 



(63) 



That coincides with the solution found in P6"| reproducing the same little rip evolution as in F(R) gravity. In this 
case the equation of state parameter is, 



\w{4>)4> 2 - v{4>) 
\u{4>)4> 2 + v^) 



2hga 



3(/i e2* E +/iie-" tE ) 2 



(64) 



Clearly this EoS parameter describes a scalar field which has a phantom behavior, where its EoS parameter < — 1 
and only tends to — 1 for i E — > 00. However, as a novelty compared with the usual phantom scalar fields (Ref. [15:]), 
here there is no finite time future singularity, but little rip evolution occurs. For this example, the time for the 
dissolution of the Sun-Earth system is given by (f5"Tj) . 

Let us now attempt to reconstruct the corresponding F(R) by applying a conformal transformation to the action 
©. It is well-known that F(R) gravity can be described by a non-propagating scalar field p"2"|) , so that the conformal 
transformation must remove the scalar kinetic term in the action (IBl: 



g^uE = tfg^v , where Q 2 



exp 



± 



d<t>y/u)(<t>) 



(65) 



The corresponding Jordan frame action yields, 



S = Id Xy—g 



2k 2 



-R 



(66) 
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This action is basically the action given in (|T2l in vacuum, which is equivalent to F(R) gravity. However, for the 
case of a phantom scalar that reproduces the little rip, it is straightforward to see that the corresponding action (prJ)) 
turns out to be complex, as the conformal transformation (|65|) is also complex due to fact that kinetic term in (|64|) 
is ne gat ive. This is the known result for the case of a phantom scalar cosmology containing future singularity (see 
Ref. [3). As in the case of phantom singular Universes, for the case of little rip, a way to reconstruct a consistent 
action in the Jordan frame can be achieved by adding an extra phantom fluid (see Ref. [17]). However, in general 
a little rip Universe described in the Einstein frame owns an inconsistent action in the Jordan frame. Thus, it is 
explicitly demonstrated that little rip evolution may be reconstructed in frames of modified gravity. 

VI. DISCUSSION 

In summary, we developed the reconstructed method for F(R) theory which admits cyclic and ekpyrotic universe 
solution. As a rule, the corresponding reconstructed F(R) action is given implicitly. Nevertheless, its expansion 
around General Relativity action with curvature corrections is always possible as is shown explicitly. The comparison 
with scalar-tensor theory which also leads to cyclic evolution is made. Our results indicate that geometrical actions 
usually considered as inflation source may lead to more complicated cyclic universe. 

We also presented the reconstruction of F(R) theory leading to little rip universe. Little rip dark energy represents 
non-singular phantom cosmology. It was proposed in Ref. [5j as viable alternative to ACDM model because it is 
consistent with observational data. Being effectively non-singular one, little rip cosmology shows the dissolution 
of bound structures similarly to Big Rip cosmology. It is demonstrated that same effect occurs for F(R) gravity 
admitting little rip solution: the corresponding time left for such dissolution is estimated. It is interesting that 
transforming little rip F(R) gravity to scalar-tensor form (the Einstein frame), the little rip universe is transformed 
to qualitatively different (Big Bang type) singular universe. The appearance of initial singularity in Einstein frame 
indicates the manifestation of dissolution time in Einstein frame. Finally, it is shown that other types of non-singular 
super-accelerating universe may be also reconstructed in F(R) gravity. Note also that the kind of solutions studied 
here can be also extended to Hofava-Lifshitz gravity (see Ref. [18|). 

Note that modified gravity theories under discussion for cyclic universe are relevant mainly for early universe 
and/or may be expanded as General Relativity action plus corrections. Hence, such theories easily pass the local 
tests. From another side, it remains to understand how natural and realistic the cyclic/ekpyrotic/little rip universes 
are. The natural recipe for this purpose is the study of perturbations. However, such investigation should be done 
for different epochs showing subsequent transition from the earlier epoch to the following oneand also with account 
of corresponding matter/radiation. This is quite complicated technical problem which lies outside of the scopes of 
current work. 
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